Mode-coupling theory of the glass transition for confined fluids 
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We present a detailed derivation of a microscopic theory for the glass transition of a liquid en- 
closed between two parallel walls relying on a mode-coupling approximation. This geometry lacks 
translational invariance perpendicular to the walls, which implies that the density profile and the 
density-density correlation function depends explicitly on the distances to the walls. We discuss 
the residual symmetry properties in slab geometry and introduce a symmetry adapted complete set 
of two-point correlation functions. Since the currents naturally split into components parallel and 
perpendicular to the walls the mathematical structure of the theory differs from the established 
mode-coupling equations in bulk. We prove that the equations for the nonergodicity parameters 
still display a covariance property similar to bulk liquids. 



I. INTRODUCTION 

Cooling or compressing a liquid usually induces a freez- 
ing transition towards a crystal, which then corresponds 
to the lowest free energy state. However, in many systems 
the phase transformation can be circumvented resulting 
in a supercooled metastablc liquid where the viscosity in- 
creases by many orders of magnitude upon mild changes 
of temperature or density. This slowing down of trans- 
port eventually leads to the glass-transition phenomenon 
where structural arrest exceeds macroscopic time scales. 
One of the grand challenges of theoretical physics is to 
provide a framework that explains the microscopic mech- 
anism and the plethora of phenomena related to the glass 
transition. 

Significant progress in this direction has been achieved 
within the mode-coupling theory of the glass transition 
(MCT), which was developed by Gotze and collabora- 
tors [J, 0| ■ The theory requires only the static structure 
as input parameter and then provides a complete descrip- 
tion of dynamic density correlations. In particular, it 
yields a strong slowing down of the structural relaxation 
upon gradual changes of the static local order, eventu- 
ally leading to a structural arrest. Thus the essence of 
the glass transition is conceived as a dynamic breaking 
of ergodicity driven by the strong correlations of the con- 
stituent particles. In the vicinity of the transition MCT 
makes a series of non-trivial predictions that character- 
ize the structural relaxation. The most prominent is the 
emergence of two scaling laws in time, a phenomenon 
that does not appear to have an analog in other fields 
of physics. The first scaling law describes the dynamics 
close to a plateau value, also referred to as nonergodicity 
parameter or glass form factor, and a factorization prop- 
erty of the space and time dependence is predicted. The 
decay from the plateau to zero obeys a second scaling law 
(time-temperature superposition principle) characterized 
by stretched relaxation functions. 

Numerous aspects of MCT have been tested success- 
fully Q; examples include depolarized- light scattering, 



which nicely displays the enhancement of a minimum in 
the first scaling regime [||^6| , colloidal glass-forming sys- 
tems @, H| exhibiting the wave-number-dependent struc- 
tural relaxation, and computer simulations on binary 
mixtures [T^j revealing scaling behavior in the vicinity 
of the plateau. 

The success of the theory for simple one-component 
systems or mixtures is encouraging to advance the mode- 
coupling approach of the glass transition also to more 
complex situations, adding new degrees of freedom, vary- 
ing the dimension or introducing confinement. The 
mode-coupling theory has been applied successfully to 
two dimensions |ll| - fl3| . Similarly, the properties of MCT 
in arbitrarily high dimensions have been discussed (l4l - 
[l7j to infer if the theory becomes of mean-field type in a 
well-defined manner. Whereas these works describe sim- 
ple liquids composed of structureless particles, already 
rigid linear [l8l - l22j and arbitrarily shaped [23[ molecules 
require the use of symmetry-adapted tensor fluctuation 
densities to account for the orientational degrees of free- 
dom. Then, the intermediate scattering function general- 
izes to a matrix-valued correlation function accompanied 
by a splitting of the currents, which introduces subtle 
new mathematical properties [24j . 

A challenge for the theory is introduced by exposing 
the glass- forming liq uid to complex geometries or ex- 
ternal potentials [25| and walls. A mode-coupling the- 
ory combining aspects of quenched disorder and interac- 
tion induced vitrification has been elaborated (26l - [29| . 
which predicts an intriguing interplay of the glass tran- 
sition driven by the strong mutual interactions of the 
fluid particles and the localization transition induced by 
the obstruction by the frozen matrix. However, there 
a divergent length scale with long-wavelength anoma- 
lies (28l . |30| emerges and a refined description requires 
concepts from critical phenomena as has been worked 
out for the Lorentz problem j3ll - l33j . 

The response of the time-dependent density correlation 
function to small local perturbations requires one to con- 
sider inhomogeneous mode-coupling equations [34| . and 
the emergence of a hidden divergent length scale has been 
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predicted. 

A great deal of experimental work and computer sim- 
ulations has been devoted to confine the liquid (HJ to a 
narrow slab to investigate the role of coopcrativity and 
dynamic heterogeneities for the slowing down of trans- 
port processes. The interaction of the liquid with the 
walls has a crucial influence on the glass transition, e.g., 
for rough walls an increase of the critical temperature in 
comparison to bulk liquids has been reported [36l - l39j . For 
smooth repulsive walls an increase, e.g., for soft sphere 
mixtures [40(, was found as well, whereas the opposite 
was identified for polymers plT - l43l ] . Experimental results 
using confocal microscopy on colloidal hard-sphere sus- 
pensions between two smooth walls reveal a decrease of 
the critical packing fraction and a slower dynamics at 
the walls in lateral direction [441-46] . Recently, an ad- 
ditional slowing down of motion has been reported due 
to an enhancement of effective surface roughness by im 
mobilized particles 



pie one-component liquids are covariant under a linear 
transformation of the time-dependent density correla- 
tor 0, l54j . This covariance, which has been proven to 
be valid for multicomponent systems as well (55j . has 
strong implications for the properties of the solutions of 
the MCT equations @, [54| [55|. The decomposition of 
the current density for confined liquids into a parallel 
and a perpendicular component leads to MCT equations 
of a different mathematical structure than for one- or 
multicomponent liquids. Here we provide a first step 
to demonstrate that some of these properties hold also 
within the mode-coupling theory for confined liquids by 
showing the covariance property of the MCT equations 
for the noncrgodicity parameters. In particular, we prove 
the existence of one of its solutions distinguished by a 
maximum principle. 



47[. Computer simulations for the H. MODEL AND INVARIANCE PROPERTIES 



diffusive dynamics of hard spheres found that the diffu- 
sivity displays pec uliar behavior as the distance of the 
plates is varied [48| , oscillating similarly to static quanti- 
ties like the excess entropy obtained by density-functional 
theory [4^]. Molecular dynamics simulations for water 
confined to silica pores have been analyzed in terms of 
the universal aspects of MCT for bulk liquids close to the 
glass transition [50l - l52j . 

Recently, wc have introduced a mode-coupling theory 
for the glass transition in slab geometry based on sym- 
metry adapted density fluctuation modes. For the case 
of a hard-sphere liquid a nontrivial interplay between the 
length scale of confinement and the average distance be- 
tween the particles has been predicted |53(. Thereby a 
facilitation of the glass transition close to half-integer val- 
ues of the distance with respect to the hard-sphere diam- 
eter was found in the theory In contrast, at distances 
near integer multiples of the particle diameter the liq- 
uid phase remains favored for higher packing fractions, 
which allows us to be interpreted as a manifestation of 
commensurability effects. 

In this article we provide a detailed derivation of the 
mode-coupling theory for confined fluids. The theory de- 
scribes the dynamics of simple fluids confined by two 
parallel flat hard walls. We discuss the residual sym- 
metries and design a symmetry-adapted complete set of 
Fourier modes to decompose the density fluctuations in 
real space. Then the density-density correlation func- 
tion is expanded in terms of a matrix-valued interme- 
diate scattering function and exact equations of motion 
are derived using the Zwanzig-Mori projection-operator 
formalism. The currents naturally split into components 
parallel and perpendicular to the surfaces. The mode- 
coupling approximation is applied for the force kernel 
leading to a set of closed equations of motion for the 
generalized intermediate scattering function. The the- 
ory requires as input the density profiles and the static 
structure factors of the confined liquid. 

The MCT equations of the glass transition for sim- 



We consider a simple liquid comprised of N identical 
particles of mass m without inner degrees of freedom en- 
closed between two flat, hard, and parallel walls, which 
are separated by a distance L. The area of the wall 
surfaces is denoted by A, and a thermodynamic limit 
A — > 00, N — > 00 is anticipated, such that the area den- 
sity no := N/A and the wall separation L remain con- 
stant. Adapted to the geometrical constraints a coordi- 
nate system is introduced such that the z axis is per- 
pendicular to the hard surfaces located at z = ±L/2. 
Here we adopt the convention that the location of the 
surfaces are taken to confine the centers of the particles 
to \z\ < L/2. In a real experiment the particles cannot 
approach the plates further than a hard core radius cr/2, 
e.g., for hard spheres the effective distance of plates is 
then H = L + a. For the development of the theory it is 
convenient to use L as the relevant confinement length. 

In the following the in-plane coordinates are abbrevi- 
ated by r — [x, y). The positions of the centers and mo- 
menta of the A^-particle system are specified by {x n } = 
{{r n ,z n )} = (xi,...,x N ) and {p n } = {(P n ,P*)} = 
(pi, ■■■,Pn) with the in-plane momenta {P n }- The po- 
sitions and momenta of the particles evolve according 
to Newton's equations of motion and the corresponding 
Hamilton function is given by 

JV -2 

H({x n }, {p n }) = £ H + V({x n }) + U({z n }). (I) 

n— 1 

For simplicity the mutual interaction between the parti- 
cles is assumed to be pairwisc additive 



N 



V({Xn}) = J2 V d f « 



(2) 



such that the two-particle interaction preserves linear 
momentum and angular momentum. The walls confine 
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the particles between the flat surfaces. Additionally, a 
specific particle-wall interaction such as adsorption in- 
duced by hydrophilic or adhesive, respectively, hydropho- 
bic or cohesive forces, can be included, 



N 



U({z n }) = '£U(z 



(3) 



where 



U{z) 



U w {z) for \z\ < L/2, 
oo for \z\ > L/2. 



(4) 



For the case of identical walls, the specific interaction 
displays the additional symmetry U w (—z) = U w (z). 

Note that the wall constraint drastically changes the 
structure and dynamics of the liquid, which cannot be 
treated within perturbation theory. Therefore, a linear 
response approach as suggested in [34| is not suitable to 
capture the induced changes. 

The physical quantities characterizing the dynamics of 
the confined liquid reflect the symmetries of the equi- 
librium distribution and the corresponding time evolu- 
tion. Whereas a bulk system is on average isotropic 
and translationally invariant, and displays space- and 
time-inversion symmetry, the walls reduce the number of 
symmetry transformations. These residual symmetries 
are determined by investigating the invariance proper- 
ties of the Hamilton function. A phase-space transfor- 
mation O is called a symmetry of the system if it leaves 
the Hamilton function invariant H({Ox n }, {Op n }) = 
H({x n y, {p n }), at least in the limit of large system sizes. 
The set of all space-time symmetries compatible with 
the constraint is generated by the elementary transfor- 
mations 



--{Tj,K z (a),V xz ,Vy Z ,{IL},T u T±}, 



(5) 



where 7j are arbitrary in-plane translations by a vector 

d = (d x ,d y ,0), lZ z (a) are rotations around the z axis by 
an angle a € [0, 2tt). The elements V xz and V yz corre- 
spond to reflections at the x-z and the y-z plane, respec- 
tively. A permutation n of the particle labels 1,2, ... ,N 
also leaves the statistical properties as well as the dy- 
namics unchanged. Furthermore, shifts 7t by a time t or 
time reversal 7± leave the dynamical laws invariant, thus 
do not change correlation functions. 

For identical walls the symmetry group is larger and is 
generated by 



^{R n (TT),l,V xv }, 



(0) 



where TZ^tt) are rotations by an angle 7r around an axis 
n = {n x , n y , 0) in the x-y plane. A space inversion is gen- 
erated by T, and V xy indicates a reflection at the dividing 
x-y plane. 



III. CORRELATION FUNCTIONS IN REAL 
AND WAVE-NUMBER SPACE 

In this section we present the correlation functions in 
real space adapted to the confined geometry. We start 
with the most basic quantity describing a liquid, the mi- 
croscopic particle density 



N 



p(f,z,t) = ^2 s [r~ r„(t)]6[z - z n (t)], 



(7) 



where the dependence on the initial position in phase 
space is omitted for simplicity. Due to translational sym- 
metry parallel to the surfaces, 7^*, the equilibrium density 
varies only in the z-direction 



n(z) = (p(r,z,t)). 



(8) 



Here the angle brackets (•) indicate canonical averaging 
over the initial conditions in phase space. We then in- 
troduce the fluctuations 6p(r, z, t) := p(r, z, t) — n(z) and 
define the density-density correlation function, which cor- 
responds to the Van Hove function [56| , 

G(\f- ?\,z, z', t) := -(5p(f, z, t)8p{r>, z' , 0)). (9) 

By translational and rotational symmetry in plane, Tj* 
and lZ z (a), the Van Hove function only depends on the 
modulus of the in-plane distance \r — f"|. Furthermore, 
by time reversal symmetry, 7± , it is an even function of 
time 



G(|r - f\,z, z', t) = G(|f - f\,z, z', -t), 



(10) 



and by time translational symmetry, 7t, it is symmetric 
with respect to interchanging the positions 



G(\r-r*\,z,z',t) = G(\r'-r\,z',z,t). 



(11) 



For the case of identical walls, the density profile is sym- 
metric, n(z) = n(—z), and the Van Hove function is in- 
variant under simultaneous reflection of z and z' 



G(\f-r'\,z,z',t) = G{\f- f |, -z, -z', t). 



(12) 



We expand all quantities in terms of symmetry- 
adapted Fourier modes. For the z direction we em- 
ploy a discrete set of Fourier modes exp(— iQ^z) with 
wave numbers = 2np/L,p 6 Z. These constitute 
a complete set exp(iQ^z) exp(— iQ^z 1 ) = LS(z — z') 
of orthogonal functions in the finite interval [—L/2, L/2]: 

/_ L ^ 2 exp(iQ M z) exp(— iQ 1/ z)dz = L5^ v . Hence the equi- 
librium density profile is expanded in discrete modes 

n ( z ) = 2^2 n ^ cxp ^' l( ^ iiZ ^ ( 13 ) 

where sums over greek subscript indices are to be taken 
over all integer numbers Z. The corresponding Fourier 



4 



coefficients are obtained as 

L/2 

n M = J dz n(z) exp(iQ A1 z). 

-L/2 



(14) 



Since n(z) is real, = n* . We shall also need the 
local specific volume v(z) := l/n(z). By the convolution 
theorem, the Fourier coefficients fulfill 

^ = *52K-k V k-» = L2 ^- ( 15 ) 

K K 

For symmetric walls the coefficients are real, 



(16) 



and by the previous relation, they are also symmetric, 
rift = n-^v^ = v-^. 

We decompose the spatial dependence parallel to the 
surfaces into ordinary plane waves, e~ lq ' r , where the 
wave vectors q = {q x ,Qy) are treated initially as discrete 
(Qx,qy) € (27r / \fA)J? . For example, the microscopic 
density is 

p(r, z,t) = j J2 \ E *) expi-iQ^e-**. (17) 

q M 

If one performs the thermodynamic limit, such that q 
becomes a continuous variable, sums are replaced by in- 
tegrals (1/A) J^q- ■ ■ -> (2tt)" 2 Jd 2 q... as usual. The fun- 
damental quantities of interest are the expansion coeffi- 
cients Pfi(q,t) called density modes, 



N 



(q, t) = J2 exp[iQ M ^(t)] e i( ^»« (18) 



with corresponding fluctuations 8p^{q,t) = p^{q,t) — 
(Pfj,(q,t)). Since {pfi{q, t)) = An^S^a the correction is 
relevant only for vanishing wave vector parallel to the 
confinement. 

Expressing the density fluctuations in real space 
p(f,z,t) by its Fourier decomposition yields an expan- 
sion of the corresponding Van Hove function 

G(\r - f|, z, z', t) =±£ -L M« t) 



A^ L 2 

q V 



x exp [i(Q M z - QuZ 1 )} e 



')] iq (f-r ) 



(19) 



in terms of an infinite matrix [S(q,i)]^„ = S liu {q,t) 1 
which generalizes the intermediate scattering function, 



S„ v {q,t) = ^(6 Pfi (q,t)*5p„(q,0)). (20) 



The translational invariance along the direction of the 
walls manifests itself in the appearance of a single wave 



vector q, whereas perpendicularly two indices arc re- 
quired. Furthermore it depends only on the magnitude 
q = \q\ due to rotational invariance in the plane, TZ z (a). 

Reversely, the generalized intermediate scattering 
function is obtained from the van Hove function by 
Fourier transform, 

L/2 L/2 

S^(q,t)= Jdz [dz' Jd{?-rt)G{\r-?\,z,z',t) 

-L/2 -L/2 A 

x exp [-i(Q„« - QuZ 1 )] e-^^'). (21) 

The initial value S flv (q) := S^(q,t = 0) characterizes 
the equilibrium structure of the fluid in the slit and will 
be referred to as generalized static structure factor. Note, 
that the hermitian matrix S(q) >z is non-negative, 
i.e., for any set of complex numbers y„, the inequality 
Em- vZSp,(q)Vv>0 holds. 

The space and time symmetries imply relations be- 
tween the matrix elements of S fJll ,(q, t). The translational 
7^-and rotational symmetry lZ z (a) have already been ex- 
ploited. By time reversal symmetry 7± and time trans- 
lation 7t, the intermediate scattering function is even in 
time and the matrix is hermitian, 

S^(q,t) = S^fa-t) = S vfl (q,t)*. (22) 

For symmetric walls, the inversion symmetry X yields 

S^(q,t) = S flv (q,t)* = S-n- v (q,t), (23) 

i.e., the matrices are real symmetric and invariant under 
simultaneous change of sign of the mode indices. 

The conservation of the particle number within the slit 
geometry is encoded in the continuity equation 



d_ 

di 



p(r,z,t) + W-j(r,z,t) = 0, 



(24) 



where the microscopic particle current density is given 
by 



JV 



j(r, z,t) = Y, ^S[r- f n (t)]8[z - z n (t)]. (25) 



By the symmetry of the system the currents split natu- 
rally into current densities parallel, 



JV 



j"(r,z,t) = J2 



Pn{t) 



8[r-f n (t)]8[z-z n {t)], (26) 



n=l 



and perpendicular to the surfaces, 

JV 



j^{r>,z,t) = Y J 



Pn(t) 



S[f - r n (t)]S[z - z n {t)]. (27) 



In this work we need only the longitudinal components 
that contribute to the particle conservation law. Then 
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the divergence V-j(r, z, t) = V?- j"(r, z, t)+V z j ± (r, z, t) 
consists of two decay channels, with currents that are 
represented in the Fourier domain, 

1 N 

n=l 

(28) 

Here we abbreviate the unit vector q = q/q and intro- 
duce the selector b a (x, z) = xS a u + z5 a ,±, which will 
simplify the subsequent manipulations. Consequently, a 
spatial Fourier expansion of Eq. (|24p leads to the conti- 
nuity equation for the density and current density modes 

d tPli {q,t) = i J2 b a (q,Q^(q;t). (29) 

a=\\,± 

It is instructive to consider also the current density cor- 
relator matrix [^T(q, t)]"P = J^{q,t), with matrix ele- 
ments defined by 

Jtf(q,t) = ^(jZ(q,ty£(q,o))- (so) 

In particular, its initial value jffi{q) = Jffijq, t = 0) 
can be evaluated explicitly (see Appendix lAj) 

^/(?) = — ^<W- (31) 
r m no 

By Eq. (|15[) its inverse matrix can be expressed in terms 
of the local specific volume 

\J-\q)\t = ^fij^^. (32) 

Applying the same reasoning as above, one easily derives 
the symmetry relations for the current-current correlator, 

J$ (ff, *) = J$ (<?> -*) = J% (q, ty, (33) 

and for symmetric walls additionally, 

J^(q,t) = J^(q,t)* = J^,(q,t). (34) 

The emergence of the channel indices (a, (3) for the cur- 
rent density correlator matrix, that represent the split- 
ting of the currents into a parallel and perpendicular 
component, occurs in the sam e sp irit as has been intro- 
duced for molecular liquids [Tj| or a single molecular 
solute [l9j]. There, the currents naturally split into a 
translational and an orientational part. 

IV. ZWANZIG-MORI 
PROJECTION-OPERATOR FORMALISM 

In this section the equations of motion for the gener- 
alized intermediate scattering function are derived with 
the hel p of the Zwanzig-Mori projection-operator formal- 
ism 0,113. 



The dynamics is driven by Newton's equations of mo- 
tion, which implies that the time evolution of phase- 
space functions A(t) = A({p n (t)},{x n (t)}) is obtained 
by d t A(t) = {A(t),H} = \£A{t), where £ is referred to 
as the Liouville operator. The formal solution then reads 
A(t) = cxp(i£t)A, where we adopt the convention that if 
no argument is provided the phase space function refers 
to the initial time t = 0. 

The set of fluctuating phase space functions is nat- 
urally equipped with a Hilbert space structure via the 
Kubo scalar product (A\B) = (5 A* SB) as correlation 
functions between fluctuations 5 A = A — (A). One eas- 
ily convinces oneself that the Liouville operator is her- 
mitian with respect to the Kubo scalar product, see, 
e.g., Q for the mathematical rigor. Dynamic corre- 
lation functions can then be represented as matrix ele- 
ments (SA(t)*SB) = (A\1Z(t)\B) of the backwards-time 
evolution operator lZ(t) = exp(— iCt). The projection- 
operator formalism relies on an exact reformulation of 
the operator identity dtH(t) = —i£1Z(t) to 

d t VK(t)V + YPCVn{t)V 

+ f dt'T£Qexp[-iQ£Q(t-t')}Q£T1l(t')T = 0, 
Jo 

(35) 

valid for any orthogonal projection operator V\ see Ap- 
pendix [B] Here Q = 1 — V denotes the projection onto 
the orthogonal complement, and 72.g(i) = exp(— \Q£Qt) 
is referred to as the reduced backwards-time evolution 
operator. 

Here we derive a formally exact equation of mo- 
tion for the generalized intermediate scattering function 
S»v(q,t) = {p^q)\K{t)\p v (q))/N. First we use the den- 
sity as distinguished variable and introduce the projector 

^^SEWS))^ 1 ^^®!, (36) 

q 

with corresponding orthogonal projection operator Q p = 
1 — V p . Sandwiching the operator identity Eq. ((35)) be- 
tween the distinguished variables, one derives the first 
equation of motion 

S^(q, t) + W K pK (q, t - t'^s-'^usxAq^ t')dt' = 0. 

(37) 

Here we observed that in Newtonian dynamics 
(Pfj,(q)\£\p v (q)) = 0, such that the second term in 
Eq. ([33)1 does not contribute. The third term can be 
simplified using Q p £\p p {q)) = £\p p (q)) and leads to the 
memory kernel K(g, t) with matrix elements, 

K pu ( q ,t) = ±(£p^(q)\ii Qp (t)\£pAq))- (38) 

In contrast to bulk systems, the current densities display 
two relaxation channels, one in plane and one in the per- 
pendicular direction. By the particle conservation law, 
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Eq. (|29|) . the memory kernel naturally splits into four 
parts, 

K^(q,t)= b a (q 7 Q p )IC a J(q,t)b^q 7 Q u ), (39) 

a/3=|l>-L 

with the reduced current-current correlation matrix 

= ^<#(S)rrc e ,(t)|i£($)>. (40) 

The symmetries of JC"@ (q, t) are identical to the current- 
current correlation function (q, t) . Its initial value 
coincides with the equilibrium static current-current cor- 
relator, Eq. pip . Solving the first equation of motion, 
Eq. (|37|) . for S lll/ (q,t) to second order in the lag time t, 
one derives the short-time expansion 

S^(Q, t) = S^(q) - iMH^ (<? 2 + Q Q v )t 2 + 0(f). 
2 m no 

(41) 

The parallel relaxation gives rise to a term for the motion 
along the plates and a second one for the flow perpendic- 
ular to the confinement. Note that different mode indices 
contribute in a nontrivial way to the decay of S pu (q, t) al- 
ready at order 0(t 2 ), due to the breaking of translational 
symmetry. 

The reduced current correlator is not suited as a start- 
ing point for approximations for the slow dynamics. 
Rather, we employ a second Zwanzig-Mori step for the re- 
duced backwards-time evolution operator 72.g p (£), hence 
we make the replacement C h-> Q P CQ P in the operator 
identity Eq. (|35[) . The new orthogonal projector is con- 
structed from the current kets \j p (q))- Since currents 
to different relaxation channels are mutually orthogonal, 
Eq. (|3ip , the projector splits into two commuting orthog- 
onal components, 

Vj = J2 V h ( 42 ) 

oHI.-L 

where the individual projections are represented by 

v t = ^EEtom^^isot®!- (43) 

q v-v 

By time-inversion symmetry the projections on the cur- 
rents are orthogonal to the densities, and the three pro- 
jection operators V p , Vj , Vj~ mutually commute. Then 
the Zwanzig-Mori procedure yields the second equation 
of motion, 

+ E E r^(«.*-*')[^- 1 (?)]ZA^(9.* , )d*' = ) 

«A ■y=\\,± J ° 

(44) 

where we observed again that S{q) is diagonal in the 
channel indices. The memory kernel 9Jt(q, t) of the fluc- 
tuating forces has matrix elements 

a>$(9,*) = ^(q)\CQe X p[-iC Q t}QL\j^q)), (45) 



where Q = Qj Q p = 1—Vj—V p projects onto the orthog- 
onal subspace spanned by the density and the currents. 
The dynamics of this subspace is generated by the re- 
duced Liouville operator Cq := QCQ. 

The exact equations of motion assume the form of 
matrix-valued integro differential equations, where the 
memory effects emerge via the convolution integrals. We 
note that due to the two decay channels both integro dif- 
ferential equations, Eqs. (|37)) and (|44j) . with first-order 
derivative in time cannot be replaced by a single integro 
differential equation with second-order time derivative, 
quite in contrast to simple one and multi-component liq- 
uids. The equations simplify in the Fourier-Laplace do- 
main, convention 

/>oo 

S ll „(q,z)=i dtS^(q,t)exp(izt), Im[z] > 0, 
Jo 

(46) 

where z constitutes a complex frequency. j58| As usual, 
one infers that the Laplace transforms are analytic func- 
tions in the upper half plane ha[z] > 0, and all singu- 
larities are concentrated on the complement From 
the definition of the Kubo scalar product, the Laplace 
transforms of dynamic correlation functions are matrix 
elements of the resolvent operator (C — z) _1 , e.g., 

S P ,(q,z) = ^{ P M\(£ - Z)~ 1 \PM), (47) 

and similarly for the other correlation functions. Trans- 
forming the first equation of motion, Eq. (|37[) . yields a 
matrix equation for S(g, z) with formal solution, 

S(q, z) = - [zS~\q) + S-\q)&(q, z)S-\q)] . (48) 

By linearity, the decomposition of K(g, t) into the differ- 
ent relaxation channels, Eq. (f3"9")l . translates directly to 
the Laplace domain, 

K„ v (q,z)= J2 b a (q,Q p )!Cf„(q,z)bP(q,Q„). (49) 
a/3=ll,-L 

Last, the second equation of motion, Eq. (|44"|). allows us 
to calculate the current kernel by matrix inversion 

it(q, z) = - [zj-\q) + J-\q)Tl(q : z)J-\q)] ~* . 

(50) 

Up to this point, all equations are exact and all features 
specific to the interactions within the liquid and the wall 
are encoded in the force kernel Wt(q,t). Close to the 
glass transition, we anticipate that forces due to inter- 
actions persist for long times, implying that the Laplace 
transform fXXt(q, z) becomes large for small frequencies. 
By Eq. (f5T)]) the current correlator K(q, z) becomes small 
in this case reflecting that transport is drastically sup- 
pressed. The first equation of motion in the Laplace 
domain, Eq. (|48| . implies that the density correlation 
function S(q, z) diverges for z — > at the glass transition 
due to the slowing down of the structural relaxation. 
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V. MODE-COUPLING THEORY 

The Zwanzig-Mori formalism expresses the density dy- 
namics in terms of the force kernel Wl(q, t). To close the 
set of dynamic equations we need to specify the force 
kernel by a suitable approximation. The basic insight is 
that caging by neighboring particles is the driving mech- 
anism for the slowing down of the dynamics. Yet, the 
caging forces entering the force kernel are generated by 
the interactions with the particles, i.e., by products of 
density modes. Here, we rely on the mode-coupling idea 
for supercooled liquids 0] to establish a connection in the 
temporal domain between the force kernel as a functional 
of the density correlation functions. The goal is thus to 
derive a microscopic theory without free parameters al- 
lowing us to evaluate the complete dynamics, including 
the long-time structural relaxation, from a set of self- 
consistent equations. 

We implement the mode-coupling idea following the 
strategy of simple bulk liquids: The forces are projected 
onto a set of fluctuating density-pair modes and the re- 
sulting four-point correlation function with reduced dy- 
namics is factorized into a product of density correlation 
functions with the original dynamics. The technical pro- 
cedure is to identify first an orthogonal projection oper- 
ator onto the pair fluctuating modes: 

V pp = \Sp(l)Sp(2))g(12;V2')(Sp(l')Sp(2')\. (51) 

11/22' 

Here we followed Ref. [l8[ to simplify the notation by 
combining the wave vectors and mode indices into su- 
per indices i = and i' = {(f il p' i ). The matrix 
5 , (12; 1'2') ensures idempotency, V 2 pp = Vpp, by the nor- 
malization condition 

£>(12; l'2')(5p(ir6p(2'rSp(l")Sp(2")) 

1'2' 

= i[5(l,l")<5(2,2") + <5(l,2")<5(2,l")]. (52) 

The essential part of the mode-coupling approximation is 
the factorization of the dynamical four-point correlation 
function into dynamical two-point correlation functions, 

(5p(iy5p(2yexp[-i£ Q t]5p(l')5p(2')) 

« N 2 [S(l,l',t)S(2,2',t) + (1' o 2')]. (53) 

Specializing to t = yields an approximate factorization 
of the static four-point correlation function 



(6p(iySp(2rSp(l')Sp(2')) 
n(p(l)\p(l'))(p(2)\p(2')) + (l'^2'). 



(54) 



For consistency, we employ the same factorization also in 
the normalization condition, Eq. (|52|) , which then allows 
us to determine 



3(12; 1'2') 



1 



AN 2 



{[S- 1 ](l,l / )[S- 1 ](2,2 / ) + (l'o2')}. 

(55) 



Collecting terms the mode-coupling procedure leads to an 
approximation for the force kernel as a bilinear functional 
of the generalized intermediate scattering function, 



[fm( q ,t)] 



1 



af) 

^ 27V 3 



x S^vi (?l , t)S^ 2U2 {q2,t)XP (q, q\q 2 ) 



(56) 

Note that, due to translational invariance in lateral di- 
rection to the walls, only wave vectors <fi and q 2 con- 
tribute which fulfill the selection rule q = q\ + q\. Here 
the complex-valued vertices X" lfla {q,q\^2) arise from 
the overlap of the fluctuating forces with the density-pair 
modes 



(57) 



The overlaps can be evaluated explicitly in terms of 
structural quantities 



{QC 3 a Jq)*5p^M)5p^m) = N 



q,qi+q2 



[b a Q- QuQM-^M + (1 o 2) 



E 



no 



■-b a (q,Q K )[s-\q)] K „s^Mqm) ; (58) 



see Appendix [C] Here, static correlations of the den- 
sity with pair modes occur, which introduces the triple 
correlation function 

Sa,^ 2 (q,qiq2) = ^(Spa(q)*Sp Pl (qi)Sp P2 (q 2 )). (59) 

In practice, the static triple correlations arc difficult 
to determine and therefore further approximations are 
introduced. In simple and molecular bulk liquids the 
convolution approximation Q has proven successful to 
describe the glassy behavior. In Appendix [D] we prove 
that the convolution approximation applied to a liquid 
confined in a slit (see Appendix [E] for details) leads to a 
similar vertex structure as found for simple and molecular 
liquids. As a result the vertices assume the compact form 



^^(Q, 01,92 



k B Tn 



TO L 2 

x \b a {t quQ^-^-^M + (1 2)]. (60) 



Here c M „ (q) are the matrix elements of the direct corre- 
lation function implicitly defined by the proper general- 
ization of the Ornstein-Zernike equation, 



no | 
L 2 



<q)], 



(61) 



with [v] pv 



[59j Direct inspection shows that 



~ P v(q) has the same symmetry properties as S pu (q). 
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The fluctuating force kernel enters the Zwanzig-Mori 
equation, Eq. (|50p . only in terms of the combination 
S~ 1 (q)fDl(q, z)J'^ 1 {q). This suggests to define an ef- 
fective force kernel A4.(q,t), which is then given by 

M$(q,t) = [J-\q)m{ q ,t)J-\ q )]f v ^ F^[S{t) 1 S{t)-q_ 

qi,q 2 =q-qi 

x S fllI/1 (qi,t)S fl2V2 (q2,t)y^ I/iV2 (q,q\q 2 y, (62) 
with new vertices 

= ^Y^ v *^ ba ^-^--^ c --^M + a 2)]. 

K 

(63) 

Note that the static inverse current correlator, Eq. (f32|) . 
is diagonal with respect to the channel index a but not 
with respect to the mode indices /i, v. The notation 
for the MCT functional J-[S(i),S(i);q] emphasizes the 
bilinearity with respect to the generalized intermediate 
scattering functions, which is a direct implication of the 
mode-coupling approximation. 

The MCT equations for confined fluids closely resemble 
the ones for molecular liquids. This motivates us to define 
a class of mode-coupling theories which is distinguished 
by multiple relaxation channels for current kernels. The 
mathematical properties proven in the next section there- 
fore hold not only for MCT of confined liquids but to all 
MCT theories belonging to this class. 

VI. NONERGODICITY PARAMETER 

In this section we introduce the nonergodicity param- 
eter, which plays a key role in glass physics for repre- 
senting the spontaneous arrest of density fluctuations. 
It is also known as the glass form factor and allows us 
to discriminate between an ergodic "liquid phase" and a 
non-ergodic "glass phase". The first subsection of this 
part provides general information about the long-time 
behavior of the self-consistent set of equations for con- 
fined liquids. In the remaining Subsections, wc prove 
certain mathematical aspects of the mode-coupling equa- 
tions. In particular, we show that an iteration scheme 
can be defined where convergence to a solution for the 
nonergodicity parameter is ensured. Furthermore this 
solution is distinguished by the property that it fulfills a 
certain maximum principle. These last subsections may 
be skipped upon the first reading of the paper. 

A. GENERAL DEFINITIONS 

In this section we show that the set of self-consistent 
equations for the intermediate scattering function can be 



solved for their respective long-time limits without solv- 
ing explicitly for the dynamics for all times. Here, we 
adopt the same approach as for bulk systems and em- 
ploy a nonvanishing long-time limit of the generalized 
intermediate scattering function 

F^(q) := lim S^[q,t) ^ 0, (64) 

t— > oo 

as the definition for a glassy state. In the current case 
they constitute an infinite matrix that inherits the her- 
mitian structure of the scattering function, 

F^{q) = F^{qy. (65) 

For symmetric walls, the inversion symmetry T yields 
additionally 

F^(q) =F^(q)* = (66) 

so the matrices are real symmetric and invariant under 
simultaneous change of sign of the mode indices. 

We argue that on general grounds the nonergodicity 
parameter is a nonnegative matrix F(q) y 0. Indeed, 
for any set of complex numbers y u , '^2 fU/ y^ l S llv (q,t)y L , 
constitutes the autocorrelation function of the variables 
y^Pfiifl> ^) an d as such its long-time limit is non- 
negative @]. 

By the Laplace transform a nonergodic contribution 
results in a zero-frequency pole, S^ v {q, z) — ~F^(q)/ z + 
(smooth), for small complex frequencies z. Reversely, the 
nonergodicity parameter can be obtained from the limit 

Fp„{q) := - lim zS^iq, z). (67) 

z— >0 

By the mode-coupling approximation, an arrest of the 
density modes is accompanied by a freezing of the force 
kernel, 

A/" (q) := M.(q, t — > oo) = J-\F, F; q], (68) 

which is again a non-negative matrix with respect to the 
double index 7 := (a, fx),6 :— {P,v) since it is the long- 
time limit of an autocorrelation function. We demon- 
strate in the next subsection that the MCT approxima- 
tion preserves this property. The simultaneous freezing of 
the forces and the densities implies again a zero- frequency 
pole for M™@(q, z) and by Eqs. (|50|) and (j49|) . the current 
correlator vanishes for small complex frequencies z — > 
as Kfiviq, z) = zG^q) + o(z) where 

GW?)= E b a {q,Q,)[M- l (q)]^(q,Q v ). (69) 

o^=|| ,-L 

In particular, one infers G(q) ^ 0, since the inverse of a 
non-negative matrix inherits the same property as well as 
its contraction with respect to the channel indices. From 
the first equation of motion, Eq. (|48p . the nonergodicity 
parameter can be evaluated as 

F(< z ) = [S- 1 ( g ) + S- 1 ( (Z )G( (Z )S- 1 ( (Z )]- 1 

=S(q)-[S-\q) + G-\q)r 1 . (70) 
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The long-time limit of the mode-coupling equations is a 
solution of the set of Eqs. (|55| - ([70| . To avoid cumber- 
some notation we allow G(q) to become formally infi- 
nite; in that case we put G _1 (<7) = 0. In general, these 
equations possess many solutions, in particular F(q) = 
represents the trivial solution, which corresponds to an 
ergodic liquid. 



B. POSITIVITY OF THE MODE-COUPLING 
FUNCTIONAL 

To demonstrate the positivity property, it is convenient 
to introduce the pair-mode indices a := (fix, fi 2 ) and b := 
iyx, V"i)- Then Eq. (|62j) allows for the compact expression 

^[F,F;g]^ 

= ^ E E y a(^ 9i&)[F(gi) ® F(q 2 )] ab y 5 b (q, M2) 

ab qi,q2=q-qi 

(71) 

where <£> denotes the Kronecker product in the space 
of mode indices, E tllUl (q 1 )F li2 ^ 2 (q 2 ) = [E(gi) (g) 

F ( < 72)]a=( Ml , AI2 ).b=(i/ 1 ^ 2 )- 

Sandwiching Eq. (|71[) between complex-valued tuples 
s 7 and summing over 7 yields 

= ^EE E ^(g, 9 lg 2 )[P( gi )®P( 92 )] a6 

75 ab 91,92=9-91 

= ^E E Zatf,m)*\F(qi)®n<b)u 

ab 91,92=9-91 

xZbfa&fa). (72) 

The first ingredient is the Kronecker product of non- 
negative matrices which is non-negative again. Second, 
a contraction with complex- valued tuples Z a (q, q\q 2 )* ■= 
s~ l *y i 2 (q, qich) is performed yielding a non-negative 
number. Thus, the mode-coupling functional maps 
non-negative matrices in the mode indices [A, v to non- 
negative matrices with respect to the double indices 7, S 
for each wave vector q. Generically, all vertices are non- 
vanishing and all components of the functional are posi- 
tive matrices ^"[F, F; q] y provided the arguments are 
positive, F(q) y 0. 



C. A CONVERGENT ITERATION SCHEME 

First, we show that the mode-coupling functional 
Af[F;q] := ^"[F,F;g] preserves the following partial or- 
dering: F y E if F(q) — F(q) y for all q. It is conve- 
nient to use a representation of Eq. (|71l) that makes the 



symmetry upon exchanging the slots manifest, 
J"[F, E; qV & = -L £ E ^ 

o-b 91,92=9—9*1 

X [F{q 1 )^F{q 2 ) + F(q 1 )(g>F(q 2 )} ab yi(q,q 1 q 2 )*. (73) 

In the following we suppress the dependence on q and 
all operations are to be understood componentwise for 
each q. For given F y and E X the arguments for 
showing the positivity of the functional are easily adapted 
to show ^[F,E] y 0. Assuming F y E one derives 
A/[F] - Af[F] = F[F + E, F - E] y 0. Thus the mode- 
coupling functional preserves ordering Af[F] y A/"[E]. 

Since inversion reverses ordering it follows that 
A/" _1 [E] — Af^F] y and therefore also the contrac- 
tions, Eq. (HH), fulfill G[E] — G[F] y 0. Eventually, the 
mapping 

* 

X[F] :=S- [S^ + G- 1 ^]]" 1 , (74) 

is continuous and also preserves the ordering 

X[F] - X[E] y 0. (75) 

Since G[F] y 0, positivity is inherited for the images of 
the mapping X[F] y 0. Furthermore S y X[F] for F y 0, 
and all fixed points F y fulfill S y F. 

We define a sequence F( n+1 ) = X[F^] with initial 
value F<°) = S y 0. Since G[S] y the first iteration 
leads to a matrix that is strictly smaller S y F' 1 '. By 
induction one infers that the sequence is monotone and 
bounded S y FW y F< n+1 ) y 0,n = 1,2... and thus 
converges to some non-negative fixed point F y 0. 



D. COVARIANCE AND MAXIMUM 
PRINCIPLE 

Here we show that the limit F obtained by iteration 
with initial condition F^ -* = S represents a maximal so- 
lution in the sense that all other non-negative solutions 
F y of the equation 

F = X[F] (76) 

arc smaller or equal than F, i.e., FXF. F is uniquely 
determined by this maximum property. The correspond- 
ing proof is based on the covariance of Eq. (f76")) under the 
linear transformation 

T:F4 T[F] = F F =: F, (77) 

which maps F to and S y to S = S - F y 0. The 
latter relation follows since Ea. (fT4")) requires S y X[F] 
for all F y 0. The requirement of covariance implies 
that there is a transformed map X[F] such that 

F=X[F]. (78) 



10 



The transformed functional X[F] is chosen such that it is 
linked to the original functional X[F] via 



X[Fl = XfFl 



(79) 



This requirement directly ensures that if F is a fixed point 
of X, then F is a fixed point of X[F]. Thus, analogously 
to Eq. ([74|) we define G[F] by 



X[F] = S - [s -1 + G _1 [I 



(80) 



Substituting X[F] from Eq. flU]) and X[F] from Eq. JgDl) 
into Eq. ((79)) and taking into account that X[F] = F we 
find the renormalized functional 



G[F] := [g _1 [F] - G _1 [] 



(81) 



For non- negative transformed functions F, i.e., for F y 
F, it follows from the previous subsection that (G _1 [F] — 
G _1 [F]) y 0. Consequently, the renormalized functional 
is again positive: G[F] y for F y 0. This in turn 
implies that all properties discussed in the previous sub- 
section remain true for X[F] on the subspace of nonneg- 
ative F, as well. Hence, the iteration of Eq. ([75| with 
initial value F^ ' = S = S — F >~ yields a fixed point 
F = limn-i.oo FW, which is non-negative F y 0. By con- 
struction of the transformed functional X[F], Eq. ([79")l . 
the diagram 



■X[F] 
r 

■XfFl 



(82) 



commutes. This property implies that the sequences gen- 
erated by the maps X[F] and X[F] are in a one-to-one 
correspondence: F" = F'"' — F. Thus, in the limit 
n — > oo the fixed points of the transformed and the origi- 
nal equation obey F = F — F y 0. Hence, the fixed point 
F is larger or equal to the assumed fixed point F: 



f y f. 



(83) 



Since the above argument applies to any fixed point 
F y 0, the fixed point F is the largest non-negative 
solution of the self-consistent equation determining the 
noncrgodicity parameter. This property will be referred 
to as maximum principle. Suppose now F* is a non- 
negative fixed point solution which fulfills the maximum 
condition Eq. (|83[) . i.e., F* y F for all fixed points of 
Eq. (f76]) . Since F is such a fixed point it is 



f* y f. 



(84) 



On the other hand we can also choose in Eq. ([531 F = F* 
since F is any of the fixed points, i.e., we obtain 



Together with Eq. (|84|) we conclude F — F*. Conse- 
quently the maximum property determines F uniquely. 

Let us note that the covariance principle used here 
is less restrictive than for simple bulk liquids. There 
it could be shown that the transformed mode-coupling 
functional is again of polynomial type 54 1. 



F y F, 



(85) 



VII. SUMMARY AND CONCLUSION 

The mode-coupling theory for liquids in confine- 
ment [HI constitutes a microscopic theory that is based 
on first principles. The scope of the theory is all two-time 
correlation functions, which can be measured experimen- 
tally by scattering methods, such as neutron, X-ray, or 
light scattering Q [5(| . The same correlation functions 
are readily obtained by tracking all particle positions as 
is performed in video microscopy [4J-[46| or in computer 
simulations (36M431 ] . Apart from being a description for 
liquids, the theory is also designed as a theory for the 
glass transition, where the structural relaxation slows 
down by many orders of magnitude. The strategy was 
first to derive a set of exact equations of motion employ- 
ing the Zwanzig-Mori formalism 0, which introduces a 
memory kernel that is a functional of all the microscopic 
details on the interaction of the particles among them- 
selves and with the walls. This functional is in general 
unknown, and the mode-coupling idea is to consider it 
as a functional local in time of the intermediate scatter- 
ing functions. The coupling coefficients are then called 
vertices and are determined from structural information 
only. 

For the case of confining parallel and flat walls, the 
fluctuating density field is expanded in a complete set 
of symmetry-adapted modes, which are continuous func- 
tions of a wave vector parallel to the planes and a discrete 
mode index for the Fourier expansion perpendicular. The 
intermediate scattering function is naturally generalized 
to a matrix-valued quantity with symmetry properties 
inherited from microscopic considerations. The breaking 
of translational symmetry perpendicular to the walls im- 
plies that the container can exchange momentum with 
a scattering probe, which is reflected in the nondiagonal 
elements of the intermediate scattering function. 

A peculiarity occurs since the currents associated with 
the density fluctuations naturally split into a compo- 
nent parallel and perpendicular to the container walls. 
This requires us to modify the structure of the equa- 
tions of motion from a single generalized harmonic oscil- 
lator to two coupled equations of motion with retarded 
friction. The same mathematical structure also occurs 
in the context of molecular liquids, where the currents 
consist of a translational and a reorientational part [l8l — 
I20I |23| . There the splitting was necessary to ensure that 
the structural relaxation dynamics is independent of the 
moment of inertia and mass of the molecule [24 1 . 

Our approach can also be employed for different types 
of confinement, such that the density modes can be ex- 
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pandcd into a complete set of geometry-adapted modes. 
In practice this can be achieved only for systems that 
display a residual symmetry such as a rectangular duct, 
a cylindrical shell or a spherical cavity, but in principle 
also for arbitrary smooth wall surfaces. 

The MCT equations for simple liquids display a se- 
ries of mathematical prop erties, which have been demon- 
strated rigorously @, [54[ . For example, the nonergodic- 
ity parameters can be obtained as the limit of a sim- 
ple iteration scheme, which is guaranteed to converge to 
a non-negative solution, as required by general proper- 
ties of autocorrelation functions. This solution is dis- 
tinguished by a maximum principle, which is obtained 
by a covariance property of the set of MCT equations 
of motion 0, [54| . For multicomponent liquids the inter- 
mediate scattering function is generalized to a matrix- 
valued quantity where the matrix indices refer to the 
different species in the liquid. The ideas of the proofs 
can then be transferred [55[ provided the notion of pos- 
itivity is generalized to hcrmitian matrices with positive 
eigenvalues. The noncrgodicity parameters of the MCT 
for confined liquids are solutions of matrix- valued equa- 
tions with a sophisticated mathematical structure due 
to the splitting of the currents. Here we have shown 
that the mode-coupling functional is positive in the ma- 
trix sense with respect to certain superindices. Then we 
suggested an iteration scheme for the noncrgodicity pa- 
rameters with initial value F'°) = S, which is monotonic 
thereby ensuring convergence. The solution F thus ob- 
tained is non-negative and fulfills a generalized maximum 
principle. The key again was to show that the struc- 
ture of the equations determining the long-time limit re- 
flect a covariance property with respect to suitable shifts 
of the nonergodicity parameters. The proofs developed 
here readily transfer also to the case of molecular liquids, 
which display the same mathematical structure. Hence, 
our paper entails important conclusions also for the well- 
established mode-coupling approach for molecular liq- 
uids [l8l - [20l [23l [60l l6l| and it is encouraging to inves- 
tigate MCT extensions for even more complex systems. 
The mathematical implications of our work demonstrate 
the robustness of the Zwanzig-Mori procedure combined 
with the first-principles MCT approach. We are confi- 
dent that also the mathematical properties of the dy- 
namic equations of motion from the bulk MCT can be 
generalized to the case of confined liquids. In particu- 
lar, one should prove the existence and uniqueness of the 
time-dependent solution and demonstrate that for over- 
damped motion the solutions correspond to pure relax- 
ations described by a superposition of decaying exponen- 
tials only. Furthermore, the long-time limit of the inter- 
mediate scattering function is expected to coincide with 
the maximal solution obtained by our iteration scheme. 
Similarly, we anticipate that all glass transition singular- 
ities in the MCT for confined systems are of the Ai type, 
specified by the classification of Arnol'd [62 1. 

Recently, striking correlations between diffusivities of 
colloidal spheres in confinement with local packing prop- 



erties have been observed and quantified in a series of 
empirical scaling properties [49j . In particular, the mo- 
bility displays oscillations as a function of the wall sepa- 
ration, which is attributed to commensurability effects of 
the packing in confinement (48j . Since the mode-coupling 
theory for confined liquids incorporates packing effects in 
terms of generalized static structure factors, it appears 
promising that our theory constitutes a microscopic basis 
for the observed empirical correlations. 

The confining walls induce strong anisotropic correla- 
tions in the liquid [63| and cannot be treated by per- 
turbation theory. In particular, the changes cannot be 
obtained as linear response to an external potential as 
has been investigated in Refs. [HI, IHj]. Our setup re- 
quires us to consider symmetry-adapted modes from the 
very beginning such that the layering and local packing 
is incorporated in suitable static quantities. 

Confinement of a liquid can also be achieved inside of a 
porous matrix where a g lass transition can occur within 
the frozen structure [26l - l29| . In addition to the slowing 
down due to caging the interaction with the disordered 
environment can lead to a localization phenomenon. In 
contrast to flat parallel walls, the disordered obstacles 
imply an additional relaxation channel in the memory 
kernel, which in mode-coupling approximation results in 
a linear coupling to the intermediate scattering function. 
Such a linear coupling is expected also for the case of 
rough walls, where corrugations open the possibility to 
exchange momentum also in the parallel direction of the 
walls. 

The mode-coupling theory for confined liquids is a mi- 
croscopic theory that does not require parameter adjust- 
ments. Hence the theory can be tested by computer 
simulations and experiments. The required input is the 
static structure which is assumed to be known. The 
MCT equations involve the three-point static correla- 
tion function, which is typically difficult to determine. 
Applying a static convolution approximation for inho- 
mogencous liquids [65| to slit geometry the vertices as- 
sume the same compact form as found for simple and 
molecular liquids. In particular, this approach reduces 
to the standard convolution approximation in the limit of 
bulk and two-dimensional liquids, respectively. While for 
three-dimensional homogeneous systems the convolution 
approximation has been proven sufficient to capture the 
key features of supercooled simple liquids, see Ref. [§6| 
for an exception, it remains a challenge for the future to 
clarify the quality of these different approaches. 

The most promising route for experimental tests are 
dense colloidal suspensions confined by glass plates [44| 
where the effects of commensurability can conveniently 
be studied. Our theory is applicable also for these over- 
damped systems provided the equations of motion are 
supplemented by friction terms accounting for the inter- 
action with the solvent and dropping the inertial terms. 
A more rigorous approach would rely on the Smolu- 
chowski operator in the first place and introduce suitable 
one-particle irreducible memory kernels as has been done 
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for bulk liquids [67j ■ These modifications affect only the 
short-time behavior; the structural relaxation encoded 
in the MCT memory kernels remains unchanged [68| . In 
particular, the phase diagram and the characteristic non- 
ergodicity parameters are identical for atomic liquids and 
colloidal suspensions. 

The MCT approach for the collective dynamics of con- 
fined liquids can be adapted to the case of tagged-particle 
motion, which is of particular interest since the self- 
dynamics is readily accessible in computer simulation and 
single-particle tracking methods on experimental sam- 
ples. In particular, the incoherent noncrgodicity param- 
eters are obtainable from a similar set of self-consistent 
matrix equations as for the collective ones, where the 
mode-coupling functional now couples to both the co- 
herent and incoherent motion [69| . Similarly, it would 
be interesting to study also the motion of the transverse 
currents and discuss the emergence of more than one vis- 
cosity due to the breaking of translational symmetry. 

Our equations allow for a direct generalization to the 
case of multicomponcnt mixtures, which is of particular 
interest, since they can be easily driven to a glassy state. 
For confinement it is even more important to suppress 
the nucleation of crystals since flat walls tend to facilitate 
the formation of ordered structures. Mixing effects fl2l . 
[7(il - [72| arise due to the presence of a new length scale 
characterizing the near order. In confinement this local 
packing competes with the layering induced by the walls 
and an even richer phenomenology is expected. 
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Appendix A: Static current density correlator 

The static current density correlation J^{q) = 
N~ 1 (j'j5{q)\ju(q)) is a diagonal matrix with respect to 
a and f3 as averages over unpaired momenta, e.g., P£P%, 
vanish. Inserting the current densities and pre-averaging 
over the momenta one obtains 

1 N 

n,m—l 

x (e if ^-^)exp(iQ, y z m )cxp(-ig Al z n )). (Al) 



Direct evaluation of the averages over the momenta yields 

(PX) = (Of" PnW- = <*nm mk B T, (A2) 

and with (pn(q,t)) = An^S-^ one obtains the explicit 
expression 

^(*) = — ^<W- (A3) 
^ m no 

Appendix B: Time-evolution operator identity 

The backwards-time evolution operator lZ(t) = 
cxp(— iCt) allows for the decomposition = 7?/p(i) + 
TZ Q (i) with K v {t) = VK(t) and K Q (t) = QTZ{t). By the 
equation of motion dtH(t) = —i£lZ(t), one obtains 

d t TZ Q {t) = -iQOl v (t) - iQCTZ Q (t), (Bl) 

which is formally solved by 

U Q (t) = c~ iQCt Q~i f e- iQC ^-^QLTZ-p(t')dt'. (B2) 

■A) 

Hence the backwards-time evolution operator can be ex- 
pressed as 

TZ(t) = VTZ(t) + c- iQCt Q- i f e- iQC ^-^Q£1Z v (l/)dt'. 

Jo 

(B3) 

The reduced backwards-time evolution operator can be 
cast in the explicitly symmetric form exp(— iQCt)Q = 
Qexp(— iQCQt)Q. Multiplying the previous equation 
from the right by V and from the left by VC one arrives 
at 

V£K(t)V = V£VK{t)V 

-i / dt l V£Qe~ iQCQ ^ t ">Q£VTZ{t')r. (B4) 
Jo 

Last, employing the equation of motion d t TZ(t) = 
—i£lZ(t), the operator identity 

d t vn(t)v + ivcm(t)v 

+ f dt'VCQe- iQCQ{t ' t ' ) QCVTZ{t')V = (B5) 
Jo 

follows, which is the starting point of the Zwanzig-Mori 
procedure. 

Appendix C: Evaluation of the overlap matrix 
element 

Here we calculate the scalar product 
(Q£j»(q)*tiPiJ.i (91)^2(92)) required for the mode- 
coupling vertex in Eq. (|58|) . With Q = 1 — Vj — V p and 
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VjlSp/j.! (<li)fiP{i2 (92)) — by time inversion symmetry, 
one obtains three contributions: 

(QCjZ(q)*6 Pfil (q-l)Sp»M)) 

= (3«(q)*[£5p^(qi)}5p»M)) + (1^2) 

- (C^{q)*V p [8p^{qi)5p^m])- (CI) 

For the first term the particle conservation law Eq. (f2T))) 
implies 

(^(g1*[^P Ml (gl)]5p M2 (9 2 )) 

= Y. bl ^Q^W^^) 5 P^)). (C2) 

Again, averaging over the momenta first, and then over 
the positions similar to Eq. (|A2|) . one obtains 

= S al S Ml+?2 ——b a (q- q\, 1)^-^(92)1^2(92)) 

fc^r - - 

= 5 ai 5 Ml+ ^N——b a {q- <fi, l)S^_ W!A12 (g 2 )- (C3) 

Here, translational invariancc implies conservation of mo- 
mentum parallel to the walls 9 = 91 + 9*2- 

Evaluating the projection on the density modes in the 
third term in Eq. (|Clj) leads to 



K,£T 

x (Spa{q)*Sp^ (qi)S Pll2 (q 2 )) 
x (5p <J {q)*6p f _ ll (q 1 )5p fJ , 2 (q 2 )), 



(C4) 



where particle conservation, Eq. (f2T))) . has been used 
again. Substituting Eq. (|31[) for the current-current 
static correlator the projected matrix element evaluates 
to 

{C:iZ{q)\V p \8pM)5 P ^M)) 



.knT 



V ^b a (q, Q K ) 

^ — <l Tin 



[S 1 ( < ?)]KrrS' (T ^ 1M2 (g,qig2) 



(C5) 



Here, we abbreviated the static three-point 
correlation function by S a . lJ , lP , 2 (q,qiq2) = 
N- 1 (5pa(q)*Sp^(qi)Sp ll2 (q2)), Eq. ®. Collecting 
terms one finds Eq. (|58|) of the main text: 



(Q£j-( 9 -)|5p Ml ( 9 l)5p M2 (g 2 )) 

=iV — 5 9,?i+?2{ &a (?-?i) Q^s^-^^M-i) + (1 <-> 2) 

~~~ E n *n-K.b a (q, QK)[S~ 1 (q)}KaS a , tllll2 (q,qiq2)}. 



n 



(C6) 



Appendix D: Vertex approximation 

In this appendix, we complete the calculation of the 
MCT vertex, using the convolution approximation in or- 
der to express the static three-point correlation function 
in terms of products of two-point correlation functions. 

The vertex after evaluating the overlap matrix ele- 
ments is given by three terms [cf. Eq. (|58|) ]: 



<%,<?i+<? 



x VK^feX^'fediP.^'fe)]^, + (1 o 2) 

x ^KMi^ei^)^- 1 ^)]^^^- 1 ^)]^^}. (Di) 

For the first two terms in the bracket, the sums over 
{p'i, p' 2 ) can be performed which leads to 

b a (q ■ 9i,Q M - M2 )[S _1 (gi)] M _ M2jMl 

+ b a (^q 2 ,Q^ pl )[S~ 1 (q 2 )]^ pi . fl2 . (D2) 

Inserting the Ornstein-Zernike equation, Eq. (|61|) . they 
can be recast to 



no 
L 2 



b a (q-q 2 ,Q 



(92) , (D3) 



where the linearity of the selector b a for a =|| and the 
selection rule q = q\ + 92 has been used. As for the third 
term, the convolution approximation (see Appendix [EJ 
gives for the triplet structure factor 

Sff^ti'^q, 9192) ~ y| E n /5i+/32+/3 3 u -/3i-A 1 S' C rA 1 (9) 

01,02,03 

Ai , A2 , A3 

x v- l 3 2 -x 2 S ( ^x 2 )^ 1 (qi)v-p 3 -\ 3 S i ^x 3 )^ 2 (q2), (D4) 

where we omit a redundant %q- 1 + ( j 2 prefactor. One can 
then successively sum out (p[, p' 2 , cr), (Ai, A2, A3), and 
k, to reduce this term to 



71 E K~ 0^~ 03^ ^Qp^+03) v ^-02^^-03- ( D5 ) 

' 02,03 



Further progress is made by making explicit the action of 
the selector and by using the linearity of Qp with respect 
to its index in the case a =_L. Eventually, performing 
the last summations over /3 2 and ^3, the third term in 
Eq. (|D1[) reduces to 



jj}b (?' QSjU— (Hi -H2 ) V fJ,-pl 1 -jj.2 



(D6) 
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and is found to cancel the first term in Eq. (|D3 
The vertex thus simplifies to 



n 



Tfl Lj 



,k B T 



n 



(D7) 

which has the same form as for simple Q and molecular 
liquids [Hill. 



related to h(ki,k2) through 

NS^(k 1 ,k 2 ) = n(k 1 +k 2 ) 

^2 E "(^ ~ Pi)"-(^2 - fii)h(pL,p2)- (E8) 



P1.P2 



Defining the local specific volume v(r) = l/n(r) and its 
Fourier transform v(k) such that 

i 2 n(fc! - - fc 2 ) = VS k k , (E9) 



Appendix E: Convolution approximation 

In this appendix based on Ref. [65[ , we report the ex- 
pression of the triplet structure factor of an inhomoge- 
neous fluid system provided by the convolution approxi- 
mation. We first discuss the general case, then specialize 
the equations to the slab geometry. 

Consider an inhomogeneous iV-particle fluid system 
enclosed in a rectangular box of volume V. Its one-body 
density and its total correlation function are denoted by 
n{r) and h(f i,?2), respectively, with the corresponding 
Fourier transforms, 

HQs) = J n{r)e itf dr, (El) 

h(k 1 ,k 2 ) = y , /i(f 1 ,f 2 )e j(Cl - Fl+4 - ?2) dr 1 df 2 . (E2) 
The triplet structure factor is defined as 

S (3) (fci, h, h) = ^(Sp(h)Sp(k 2 )Sp(k 3 )), (E3) 



with 



and 



6p(k) = p(k) - (p(k)) = p(k) - h(k) (E4) 



N 



(E5) 



p(k) = Y^e ik ^, 

3 = 1 

where Xj is the position of the jth particle. Note that, at 
variance with the main text, the definitions of the struc- 
ture factors in this appendix do not involve any com- 
plex conjugation to preserve the symmetry of the working 
equations. 

Following Rajan et al. (65j . the convolution approxi- 
mation for S^ 3 '(ki, k 2 , £3) reads 

S( 3 >(£i,£ 2 ,fc 3 ) « — g h & +P2+P3) 



K X ,K2,K 3 

Pl,P2,P3 



xY[n(ki + Ki) [h(-Pi, ~Ki) - h(h -pi,-ki - K t ) 

i=l 

(E6) 

This result is most conveniently reformulated in terms of 
the pair structure factor 

^ (2) (fci,fc 2 ) = ^(Mfci)^(fc2)), (E7) 



Eq. (|E8|) is easily inverted to yield 

h(ki,k 2 ) = -v(ki + k 2 ) 

M ^2 E 5 & - pi)v(k 2 - p 2 )S^(p 1 , P 7 2 ), (E10) 



V 2 



Pl,P2 



which can be injected into Eq. (|E6|) . The summations 
over K\, K 2 , and K3, can then be explicitly performed 
and, using the fact that for a closed system (e.g. canoni- 
cal system) S^ 2 \k, 0) = for any k, it follows that in the 
convolution approximation 



S< 3 >(fei,&,&) 



Y HP1+P2+P3) 



Pl,P2,P3 
Il,l2,l3 

3 



Jv{-Pi-k)S (2) (Ti,ki). (Ell) 



One can readily check that this expression reproduces 
the standard result for bulk systems. Indeed, one then 
has h(k) = N6 U;S , v(k) = {V 2 /N)5^ and S^(p,k) = 

s ( k ) s p+%,& so that 

S^(k u k 2 ,k 3 ) « <5jE 1+fe+j6l>ff 5(fci)S(fe a )S(A 3 ). (E12) 

Application to the slab geometry is just as straight- 
forward. One simply has to set V = LA, split 
each sum over a wave vector k = (q, Q^) into one 
over a transverse index p and one over an in-plane 
wave vector q, and replace h(k), v(k), S^(ki,k 2 ), and 
S^(ki,k 2 ,ka), with n^AS^, v^AS^, Sjff^ (qi)S ?1+?2 g, 
and Sltinzus^qi, q 2 , qs), respectively Eventually, one gets 



<?l+<?2+<?3,0 £6 

3 



e "ft+fe+fe n^-^-^^i.fe)- (E13) 



Pl,P2,P3 

Ai , A2 , A3 
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